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mass  of  cold  e3q>losive*  The  time  to  ignition  of  the  esqplosive  is  detexudned 
by  numerioal  methods,  and  the  way  in  idiich  this  time  varies  with  the 
px«pex*ties  of  the  hot  iqpot  and  ea^losive  is  investigated*  Four  speoiflo 
hot  spot  materials  are  considered:  copper,  steel,  glass  and  a  further 
material  with  thermal  properties  similar  to  those  of  the  explosive*  The 
first  three  of  these  were  surrounded  by  RDX,  but,  if  certain  oonditioos 
are  satisfied,  the  results  are  explicable  to  other  e^qplosives* 
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1 .  lilTRODUCTION 


In  this  report,  the  initiation  of  an  explosion  hy  a  small  hot  sphere 
of  chemically  inert  material  embedded  in  a  mass  of  cold  explosive  is 
considered.  The  object  of  the  present  calculations  is  twofold;  first, 
to  set  up  and  programme  the  appropriate  equations  in  a  finite  difference 
^orm  suitable  for  numerical  integration  on  the  high  speed  computer  AMOS; 
and  second,  to  evaluate  the  time  to  ignition  for  a  given  hot  spot  and 
explosive  and  to  investigate  the  effect  of  the  thermal  properties  of  the 
hot  spot  material  on  the  Ignition  time.  In  order  to  cover  a  wide  range  of 
explosives  an  enormous  amount  of  con?)utation  woxild  be  required;  this  did 
not  seem  justified  at  present  in  view;  of  several  approximations  and 
assvmptions  involved  in  the  mathematical  model.  Only  one  explosive,  RDX, 
has  therefore  been  treated  in  detail  but  it  should  be  pointed  out  that 
the  programnes  constrtKsted  for  AMOS  are  quite  general  and  any  fiarther 
oases  of  particular  interest  could  be  easily  examined. 

A  series  of  experimental. investigations  of  hot  spot  initiation  has 
been  carried  out  by  Bowden  and  his  co-workers  amd  details  are  given  in 
Bowden  and  Yoffe  [l].  Their  general  conclvisions  concerning  the  initiation 
of  explosion  by  grit  particles  were  that  a  particle  of  about  I0"^cms. 
diameter  at  about  500°  C  was  capable  of  causing  initiation  in  a  time  of  the 
order  of  10"'*  secs,  in  a  small  saii5)le  of  explosive.  Roughly  sx>eaking, 
the  same  orders  of  magnitude  have  been  found  in  the  present  calculations 
although*  as  we  shall  show  later,  the  effect  of  the  thermal  properties  of 
the  hot  spot  is  quite  pronovoiced.  As  far  as  the  authors  are  aware,  no 
previous  theoretical  attempt  to  study  the  effect  of  a  solid  inert  hot  spot 
has  been  made.  Rideal  and  Robertson  [2]  have  considered  initiation  by  an 
explosive  hot  spot  and  have  evaluated  critical  radii  and  tenperatures 
for  several  explosives;  some  further  work  on  explosive  hot  spots  was 
carried  out  on  the  Aiken  computer  at  Harvard  and  is  reported  by  Seager 
[3].  An  approximate  discussion  of  initiation  in  a  plastic  explosive  due 
to  the  rapid  compression  of  enclosed  air  bubbles  has  been  given  by  Coombs 
and  Thorniu.il  [4]. 

The  present  report  is  the  third  of  a  series  in  thermal  explosion 
theory.  In  the  first  of  these  L5l»  one  face  of  a  slab  of  cold  explosive 
was  imagined  to  be  placed  in  contact  with  a  constant  temperature  bath. 

The  time  to  ignition  was  calculated,  and  from  the  results,  enpirical 
relations  connecting  the  ignition  time  v/ith  the  initial  uniform  tenpera- 
ttire  of  the  sample  and  the  bath  temperature,  were  derived.  In  the 
second  [6],  a  pulse  of  high  intensity  thermal  radiation,  such  as  that 
delivered  from  an  atomic  weapon,  was  assumed  to  fall  on  a  steel  plate 
covering  a  mass  of  explosive;  critical  conditions  under  which  explosion 
only  just  uanages  to  occttr  were  evaluated  for  RDX.  The  mathematical 
details  gathered  in  paras.  3  and  4  may  be  passed  over,  if  desired, 
without  loss  of  continuity. 

2.  P0RI.1ULATI0N  OP  THE  PROBLEIjI 

Consider  a  small  hot  sphere  of  radius  a  and  initial  uniform 
tenperatxjre  Tq  embedded  in  an  infinite  mass  of  cold  explosive  at 
initial  xmiforra  tempera ttare  Ti.  This  model  is  suitable  for  quite  small 
samples  of  explosive  since  the  explosive  is  a  poor  conductor  arid,  in 
the  short  times  involved,  the  heat  penetration  is  very  small.  On  the 
assumption  that  the  thermal  properties  of  the  hot  spot  material  «uid  the 
explosive  do  not  vary  with  temperature  and  that  the  decomposition  of  the 
explosive  can  be  represented  by  a  single  chemical  reaction  of  the  zero 
onder,  the  governing  equations  can  be  vrritten. 
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Tuhere  svif fixes  *  1 '  and  '  2'  refer  to  the  hot  spot  material  and  explosive 
respectively.  In  (1)  and  (2)  the  following  notation  is  employed: 

k  =  thermal  conductivity, 

P  =  density, 

c  s  specific  heat, 

T  =  temperature , 

E  =  activation  energy  of  the  reaction, 

Q  s  heat  liberated  by  chemical  reaction  per  gram, 
and  Z  =  frequency  factor. 

Equations  (I)  and  (2)  have  to  be  solved  under  the  foUovidng  initial  and 
boundary  conditions: 


3Ta 


t  s  0;  0<  r  <  a; 

r  >  a; 
t  >  0;  r  =  a; 


r  00  ; 
r  =  0; 


Ti  =  To  , 
Tz  =  Ti  , 
Ti  a  Tz . 

dTi  _  .  3T2 

-§5:  -  kz  , 

T2  -  Ti, 

Ti  finite. 


(5) 

(4) 

(5) 

U) 

(7) 

(8) 


It  is  convenient  to  introduce  the  dimensionless  variables  t  and  0 
defined  by 
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-  —  T  F  —  /RQZPz  \/4  ^fj-7  _  _  /RQZ\  ^q-i4. 

E  1cr  r,  r  =  (— )  1(r  t. 


(9) 


E  '  kzE  '  •  -  'C2E'' 

=  aiT,  =  azr,  =  ast, 

where  we  suppose  ^  has  the  value  a  when  r  =  a. 

The  resulting  equations  can  be  ftcrther  simplified  by  writing: 

u  =  €0, 

and  it  is  also  convenient  to  introduce  the  variable  S,  defined  in  the 
interval  0  <  ^  <  a  by-  the  equation 


^  J  K? 

where  K  is  the  thermal  diffusivity  k/pc,  and  S  =  ^  when  € 
Equations  (1)  to  (8)  then  become: 
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r  >  0;  €  a  a* 

U1  a  Ua , 
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€  *  0; 

Ui  a  0* 

(80 

It  has  been  stated  above  that  the  thermal  properties  of  both 
materials  have  been  assxmed  to  be  Independent  of  temperature*  This 
asstimption  is  mathematically  convenient  but  it  should  be  pointed  out 
that*  in  general,  data  on  the  variation  of  these  properties  with 
temperature  is  rather  soaroe*  It  has  been  further  assumed  that,  over 
a  vide  tenperature  range,  the  decomposition  can  be  represented  by  a 
single  chemical  reaction  of  the  zero  order*  While  this  is  undoubtedly 
an  oversimplification,  it  is  probably  sufficiently  realistic  to  enable 
at  least  the  qualitative  features  of  the  behavioiir  of  the  system  to  be 
predicted.  It  is  further  assumed  that  convection  effects  arising  from 
the  liquefaction  of  the  explosive  can  be  ignored  and  that,  at  all  times, 
there  is  no  contact  resistance  at  the  e3Q>losive/hot  spot  interface* 

3*  FINITE  DIFFERENCE  METHOD 

Providing  a  certain  restriction  is  observed,  a  stable  finite 
difference  form  of  equations  (1  ')  and  (2^)  emerges  if  forward  differences 
are  used  for  the  time  derivatives,  and  central  differences  for  the  space 
derivatives*  With  the  definitions: 


a  ui  (nAS,  JAr), 
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and  So 
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equations  (1^) 

and  (2^)  become: 
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where 


=  O’l  (unl},d  +  -  2ui’J  )  +  0  <  n  <  IT  (10) 


Un  ,  J*i 

=  <ri 

(iJnl 

(0 

V4/>1 

=  0*2 

(ui* 

10^*  (a  +  m&O  At  exp  f-^a  +  inA€)/tim*^J»  m  >  0  (11) 


At 


0-,  =  “V(AiS)S  0-2  =  ^7(A€)*  . 


In  the  numerical  work  A^  and  AS  were  al'.vays  given  identical  values,  so 
that 

ffi  s  0*2  s  c. 


There  is  still  a  choice  of  methods  that  may  he  used  to  expand  the 
boundary  condition  (6')  in  finite  difference  form*  Here  two  methods  are 
developed,  and  compared  enipirically.  In  the  first  both  terms  are  expanded 
by  Taylor*  s  theorem  and  central  differences  used  on  the  resulting 
expression.  This  method  yields  the  boundary  equation: 


UN,  j  =  uo  ,  j  =  Ui  x^i  +  Ua  Ui  ,  j  -  Ai  UN-2,  j  -  Aa  Ua  ,  j. 


(12) 
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Alternatively  the  "hypothetical  point"  method  may  be  used  in  which  it  is 
assumed  that  each  material  is  extended  over  the  actual  interface  by  an 

(0* 

amount  A€,  the  value  of  u  at  these  points  being  defined  as  ,  j  in  the 

(O* 

case  of  the  hot  spot  material  and  u.i.,j  for  the  explosive.  The  equations 
governing  the  behaviour  of  the  hypothetical  points,  corresponding  to 
(10»  (2')»  (50  and  (6'),  are 
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UN,>i  (uN^i  ,;5  +  uiJ-1  -  2uN,j)  <r  +  UN,j, 


(0 

Uo,.>i 


(a)  (a)*  (a)  (a)  ^  (a) 

=  (ui,j+  11.1  ,j  -  2uo,j)  o-  +  uo,j  +  oAr  10^^e^Uo,j  , 
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2  N  (O  f  (O*  (*)  V 

-  (A0(1  -  v)  H  (^♦i  ^  j  -  Uff., 
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Eliminating  u*  yields  the  boundazy  condition 
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where 
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The  remaining  boundary  equations  in  finite  difference  form,  to  be  used 
with  either  (12)  or  (l2"),are 


(a) 

^,6  (“  +  n>d€)  ®i> 

vO  - 

a  0, 

iL*  ^  m  vt/t^  ®® 

Un  0  a  IMt  ~  > 
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and 


m  eo 


0  <  n  <  N 


m  >  0 


(13) 

(1*) 

(15) 

(16) 


Once  an  initial  line  is  supplied, equations  (10)  -  (14)  are  stifficient, 
for  the  evaluation  of  u  at  all  later  times.  However,  equations  (15) 
and  (16)  give  a  discontinuity  on  the  initial  line  at  the  interface  which 
is  \indesirable  at  the  commencement  of  the  computation.  To  overcome  this 
difficulty  a  formal  solution  of  the  conduction  equations  in  the  plane  case, 
and  without  the  exponential  term,  may  be  obtained  in  the  form; 


Un 


(')  .  nllMi 


»•  -  TTTS  ^  1 

2  Vj?  J 


(17) 


um,;J  =  (a  +  mA?)  I  0i  +  -■  ^  (0©  -  0i)  erfc  — 2_-  T 

L  “  ■  2/3?  J 


In  general»one  is  justified  in  starting  the  computation  from  a  small  value 
of  j  using  equations  ( 1 ?) ,  since  initially  the  temperature  in  the  explosive 
is  sxrfficiently  low  for  the  rate  of  chemical  reaction  to  be  negligible 
and,  for  a  sufficiently  large  hot  spot,  there  is  little  difference 
between  the  solutions  of  the  plane  and  spherical  cases. 
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4.  MATHEf^ATICAL  DISCUSSION 


Typical  temperature  profiles  based  on  the  equations  (10)-(17)  are 
shown  in  Pig.  1»  Ignition  was  said  to  have  occtirred,  and  the  time  t 
to  ignition  derived,  when  the  rate  of  increase  of  temperatvire  at  any® 
point  in  the  explosive  accelerated  rapidly. 

Before  investigating  the  effect  on  t  of  the  physical  parameters  a. 
To  >  ,  M  and  v,  it  is  of  interest  to  determine  the  way  in  which  these 

results  will  depend  on  the  choice  of  the  mathematical  parameters  At, 
and  o’ ,  and  also  on  the  choice  of  the  boimdary  equation. 

As  it  is  not  possible  to  compare  directly  the  numerical  solutions 
with  the  formal  solution  of  the  general  equations, we  consider  the  simpler 
case  of  a  spherical  hot  spot  embedded  in  a  non-explosive  material  and 
with  the  condition  /i  =  v  =  1,  The  solution  in  a  suitable  form  is 
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=  00 


nA€  -  t  (00  -  0i) 


[■ 
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When  the  numerical  solutions  were  compared  with  those  given  by  equations 
(18)  and  (19)  it  was  found  that  nowhere  did  the  ei*ror  exceed  0,1%  with 
the  use  of  either  boundary  eqtiation  and  that,  in  general,  the  actual 
error  was  approximately  halved  when  (12^)  was  used.  It  will  be  seen 
from  Pig.  2  ,  where  Tg  is  plotted  against  At  for  both  boundary  equations, 
that  equation  (12)  led  to  faster  ignition  ifidiile  both  equations  give  a 
common  value  of  as  Ar -*■  0,  Thus  the  use  of  equation  (12')  as  the 
boundaury  condition  saves  much  computation  time  while  the  dependence 
of  Te  on  Ar  appears  to  be  linear, making  extrapolation  to  the  exact 
solution  simpler.  The  results  quoted  later  are  derived  by  use  of  equation 
(12')  in  preference  to  equation  (.12), 


It  is  also  apparent  from  Fig.  2  that  the  ignition  time  is  sensitive 
to  changes  in  Ar.  However,  as  demonstrated  in  Fig.  3  where  results  are 
obtained  for  two  values  of  Ar,  the  form  of  the  cvirves  remains  the  same. 
Since  this  is  expected  to  be  generally  the  case,  and  in  this  report  we 
are  more  concerixed  with  the  behaviour  of  t^  rather  than  exact  numerical 
results,  only  one  value  of  Ar  has  been  taken  for  any  derivation  of  t  . 

In  any  case,  there  is  some  considerable  doubt  as  to  the  accuracy  of  ?he 
published  figtires  used  for  obtaining  the  parameters  ai ,  a* ,  as  of 
equation  (9).  If  accurate  ignition  times  are  required  these  may  be 
obtained  by  determining  te  for  two  or  three  values  of  Ar  and  extrapolating 

to  zero  in  the  manner  of  Pig.’  2. 

.-.-r  * 
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It  has  been  shown  by  O’Brien,  Hyman,  and  Kaplan  [7]  that  linear 
equations  put  into  finite  difference  form  in  the  way  described  in  para* 
3.  are  stable  if  o’  <  %  and  it  is  known  that  the  truncation  error  is  at 
a  minimum  if  O’  =  .  For  craaparison  two  values  of  o’  were  investi^ted, 

O’  =  /a  and  o’  =  %  while  At  was  kept  constant.  The  resulting  valxies  of 
te  differed  by  less  than  0*5% »  and  since  it  is  convenient  when  consid¬ 
ering  small  radii  and  relatively  large  times  to  choose  o’  as  large  as 
possible,  the  value  o’  =  %  was  used  in  all  fxjrther  calculations. 

3.  RESULTS 

Owing  to  the  large  number  of  parameters  to  be  considered,  machine 
coii?)utation  was  confined  to  give  only  specimen  results,  representati’ve 
of  a  wide  field.  For  this  purpose  four  actual  hot  spot  materials  were 
selected  with  thermal  properties  ranging  from  those  similar  to  the 
e:q)losive  itself  to  those  of  copper.  The  remaining  computation  time 
was  devoted  to  studying  the  general  behaviovir  of  the  time  to  ignition 
te  with  varying  radius  of  hot  spot  a,  initial  tenperatxires  Tq  and  Ti 
of  the  hot  spot  and  e:qolosive  respectively,  and  the  thermal  parameters 


The  variation  of  t^  with  either  T^  or  T^  is  similar  in  appearance  to 

that  of  t  with  the  radius  a.  The  gradient  for  small  initial  temperatures 
is,  hovTevSr,  not  so  steep,  and  in  the  ranges  considered  the  graphs  suggest 
a  relation  of  the  form 

log  t  =  ^T'  +  Bi  ,  (20) 

6 

where  A  and  B  are  constants  and  T'  is  either  or  T^.  This  relation 
has  also  been  found,  though  under  different  circvonstances,  by  Henkin  [8] 
and  Cook  [5]» 

Since  the  thermal  parameters  ki  ,  ka ,  Ki  ,  Ka  appear  only  in  the  ratios 
g  ,  V  it  is  sxtfficient  to  investigate  the  dependence  of  t  on  these  two 
variables.  It  is  apparent  from  Fig,  4  that,  for  small  valuSs,  r  has  an 
almost  equal  dependence  on  either  v  or  g-’  ,  or,  in  terms  of  thS  hot  spot 
properties,  on  either  Ki"^  or  ki  .  For  large  values  the  dependence  on 
ji"*  d.omimtes.  Another  fact  that  emerges  fran  Fig.  4  is  that  if 
accurate  ignition  times  are  required  for  a  hot  spot  'with  thermal 
properties  similar  to  those  of  the  explosive  (i.e,  m~1,  v-'I),  more 
accurate  knowledge  of  the  explosive  properties  is  needed 'than  if  the 
hot  spot  were  made  of  a  material  such  as  steel. 

The  dominant  variable  is  shown  more  clearly  in  Fig.  5  where  r 
is  plotted  as  a  ftmction  of  m  under  the  condition  g  =  w.  For  small 
values  of  M,  even  though  the  hot  spot  conductivity  is  extremely  large, 

T  approaches  the  adiabatic  ignition  time  (i.e,  the  time  taken  for  a 
itrge  mass  of  explosive  to  ignite  spontaneously).  Here  a  relation  of 
the  form  (20)  exists  where  M  is  substituted  for  T',  At  the  other  end 
of  the  scale  it  is  apparent  that? if  ki  approximates  to  zero, it  is  still 
possible  to  obtain  a  small  ignition  time  by  supplying  sufficient  initial 
heat  * 

The  four  representative  materials  chosen  for  the  hot  spot  were 
copper,  steel,  glass  and  a  further  material,  substance  A,  having  thermal 
properties  identical  to  those  of  the  explosive.  The  thermal  properties 
assumed  for  the  first  three  acre  listed  in  Table  1 . 
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Table  1 

Material 

ki  ( cal,  cm”^,  sec"J  (^C)"^) 

Ki  (  on?,  seal  ^  ) 

Copper 

0.85 

0*81 

Steel 

O.O85 

0*081 

Glass 

0.0025 

0.006 

The  conductivity  and  ^ffusivity  of  the  explosive  were  taken  to  he 
7  X  10“^cal.ciisecTn  and  10"*  cm.*  seoV  respectively  in  sQ.!  cases* 

Figs,  6  and  7  show  typical  ciarves  for  ignition  times  against  initial  hot 
spot  temperatures,  the  radius  of  the  hot  spot  being  O.O38  cms.,  the 
initial  explosive  temperature  50°C  and  the  e:q)losive  RDX,  The  curves 
again  approximate  to  a  relation  of  the  form  of  equation  (20), 

Since  ti  and  v  are  the  only  thermal  parameters  to  appear  in  the 
dimensionless  formulation  of  the  problem  it  follows  that, if  these  may  be 
assumed  constant  for  different  e:q>losives, results  applicable  to  these 
e:q>losives  may  be  derived  Toy  choosing  ai ,  ag,  aj  of  equation  (9) 
accordingly*  The  valvies  of  these  parameters  adopted  in  this  report 
for  three  explosives  which  approximately  satisfy  this  condition  are 
(Hideal  and  Robertson  [2]): 

Table  2 


RDX 

FETN 

HMX 

10^  a, 

0*418 

0.425 

0*377 

1700 

4700 

8983 

aj 

2770 

20800 

53800 

Fig*  8  shows  the  dependence  of  t  on  the  initial  hot  spot  temperature  To 
for  each  of  the  three  explosives.  For  the  c^arves  labelled  I,  II  and  III, 
appropriate  data  are  listed  in  the  following  table. 


Table  3 


Graph 

Explosive 

• 

Radius  of 

Hot  Spot  (oms."' 

Initial  temp, 
of  Explosive  (°C) 

Time  scale 

factor  (s) 

I 

RDX 

0*5 

15 

10* 

II 

PETN 

0*181 

11.6 

10* 

III 

HMX 

0.095 

i»6*3 

10* 

The  graphs  in  which  t  is  plotted  as  a  function  of  To,  T.  or  a,  all 
have  a  similar  appearance  In  that, after  a  certain  point  is  reached, a  aman. 
decrease  in  one -of  the  independent  variables  yields  a  large  increase  in 
t  .  This  leads  one  to  ask  if  there  is  not  a  small  inteznral  over  which  it 
might  be  said  that  at  one  end  ignition  is  due  to  the  hot  spot  at  the 
other  to  a  uniform  rise  in  the  e:^losive  ten^ieratiare .  In  this  manner  it 
would  be  possible  to  define  a  "critical"  value  of  T  ,  T.  or  a  for  a  given 
solid  hot  spot  similar  to  that  defined  by  Rideal  anS  Rotertson  [2]  for 
an  explosive  hot  spot.  They  define  "critical  temperature"  to  be  that 
initial  hot  spot  temperatiire  at  which  heat  evolution  due  to  chemical 


c  s  ) 


RRSTRICTED 


RBSTRICTED 


reaction  within  the  hot  spot  in  a  time  t,  just  balances  the  heat  lost 
by  conduction  to  the  surrounding  explosive  in  the  same  time.  This 
condition  is  given  by  the  relation 


%  «a*pt  QZe  '  ^RTo  = 


kn"*  Ta  pcdr» 


a 


the  temperatvire  Ta  on  the  ri^t  hand  side  being  determined  from  the 
linearised  heat  conduction  eqmtion  on  the  assumption  that  heat  evolution 
in  the  cold  e:qplosive  is  negligible  dxnrLng  the  time  t.  Thus  for  a  given 
value  of  t»  a  '’critical"  value  of  T^  may  be  derived.  In  the  present  ease 
however  such  a  definition  is  impracticable  and  a  graphical  investigation 
appears  to  be  necessary.  So  far  as  the  coii;>utation8  have  been  conducted 
the  curves  do  not  differ  significantly  from  equation  (20),  With  the 
present  formulation  of  the  problem, cm  investi^tion  of  the  higher  parts 
of  the  ctorves  would  involve  an  enormotis  amount  of  conqputation  and  it 
has  notf  at  present, been  attempted. 
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FIG  I.  TEMPERATURE  PROFILES-  SUBSTANCE  A  6«0-038  CMS.) 
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FIG.  2 


FIG. 2  IGNITION  TIME  AS  A  FUNCTION  OF  INTERVAL 


RESTRICTED 


RESTRICTED 


FI6.3 


FIG.  3  IGNITION  TIMES  FOR  VARIOUS  RADII 
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FIG. 4 


FIG. 4  IGNITION  TIMES  FOR  VARYING  RATIOS  OF  THE  CONDUCTIVITIES. » 

AND  DIFFUSIVITIES.it*.  OF  THE  EXPLOSIVE  AND  HOT  SPOT 
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FIG.5 


FIG. 5  IGNITION  TIMES  WITH  THE  RATIO  OF  THE  CONDUCTIVITIES  OF  THE  HOT 

SPOT  AND  EXPLOSIVE  EQUAL  TO  THAT  OF  THE  DIFFUSIVITIES  (i.t. /izv) 
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FIG.  6  IGNITION  TIMES  OF  RDX  FOR  TWO  TYPES  OF  HOT  SPOT 

RESTRICTED 


5. 10  SECS. 


RESTRICTED 


FI6.7 


FIG. 7  IGNITION  TIMES  OF  RDX  FOR  TWO  TYPES  OF  HOT  SPOT 
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FIG.8 


FIG.S^IGNITION  TIMES  FOR  THREE  EXPLOSIVES  WITH  A  STEEL  HOT  SPOT 
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